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^ : Abstract 

o> 

X I A measurement of inclusive semileptonic branching fractions of b hadrons produced in Z 

I ■ decays is presented. An enriched Z° -^ bb sample is obtained with a lifetime flavour-tagging 

^' technique. The leptonic events are then selected from this sample, and classified according to 

^ , their origin, which is determined by comparing the distribution of several kinematic variables 

using artificial neural network techniques. Using 3.6 million multihadronic events collected 

k>( ; with the OPAL detector at energies near the Z^ resonance, the values 

^ : BR(b -> iX) = (10.83 ± 0.10 (stat.) ± 0.20 (syst.) t^l (model)) % 

BR(b ^ c -> £X) = (8.40 ± 0.16 (stat.) ± 0.21 (syst.) toU (model)) % 

are measured, where b denotes all weakly decaying b hadrons and i represents either e or /i. 
The second error includes all experimental systematic uncertainties whereas the last error is 
due to uncertainties in modelling of the lepton momentum spectrum in semileptonic decays 
and b quark fragmentation. The average fraction of the beam energy carried by the weakly 
decaying b hadron, (xe), is measured to be 

(xe) = 0.709 ± 0.003 (stat.) ± 0.003 (syst.) ± 0.013 (model) 

where the modelling error is dominated by the choice of b fragmentation model. The agree- 
ment between data and various semileptonic decay models and fragmentation functions is 
also investigated. 
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1 Introduction 



Measurements of the seinileptonic branching fraction of b hadrons are important in testing 
our understanding of the dynamics of heavy quark physics and are also important inputs for 
other b physics analyses. These measurements can also be used to extract the CKM matrix 
element Vcb- Recent QCD calculations which include higher-order perturbative corrections 
|jl|, have lowered the predicted value of the semileptonic branching ratio for B mesons, 
BR^^, and now adequately reproduce the experimental results 0. These calculations also 
predict a value for (nc), the average number of charmed hadrons produced per B meson 
decay, which is consistent with experimental data 0]. 

While theoretical calculations are now in better agreement with experiment, the measure- 
ments of BR5J;, obtained at the T(4S) and Z° resonances have slightly disagreed for some 
time. The semileptonic branching fraction for B mesons has been measured at the T(4S) 
resonance to be BR^^ = (10.45 ± 0.21) % |^, whereas a combination of all results obtained 
at the Z° resonance gives BR^^ = (10.99 ± 0.23)% 0, where the superscript b indicates 
that the high energy data correspond to a mixture of B^, B°, Bs and b baryons, as opposed 
to B''^ and B° only at the T(4S) resonance. Assuming the semileptonic width, Fgi, to be the 
same for all b hadrons, as suggested by the result of [Q, and given that the semileptonic 
branching ratio is related to the lifetimes, r, by BR^^ = Fsi/Ftotai = ''"rsi, 0116 obtains 

BR5L = BRg^. (1) 

To remove the difference between the T(4S) and 7P results would require that t-q/t\^ be 0.948, 
whereas current lifetime measurements yield t-q/t^o greater than one 0. After applying the 
correction factor Te/rb, there is a difference of about two standard deviations between the 
T(4S) and 7P results for the branching fractions. 

This paper describes the measurement of the semileptonic branching fraction for all b hadrons 
produced at or near the 7P resonance using identified electrons and muons in an enriched 
Z° -^ bb sample. The measurement differs from previously published measurements 0- 
H in the use of a method designed specifically to achieve a precise determination of the 
semileptonic branching ratios rather than extracting them from a multi-dimensional fit for 
several electroweak parameters. The dependence on the semileptonic decay models has been 
substantially reduced and the correlation with the -Rb measurement eliminated. The analysis 
also allows the determination of BR(b ^ c ^ £X) and (x^), the mean fraction of beam 
energy carried by the weakly decaying b hadron, both of which are important inputs needed 
for other heavy flavour measurements, such as those described in . The agreement between 
data and various semileptonic decay models and fragmentation functions is investigated in 
the Appendix. 



2 The OPAL detector, data and Monte Carlo samples 



The OPAL detector is described in reference |10]. The central tracking system is composed 



of a silicon microvertex detector, a precision vertex drift chamber, and a large-volume jet 

1 



chamber surrounded by a set of drift chambers that measure the 2;-coordinate.[] Charged 
particles are identified by their specific energy loss, dE/dx, in the jet chamber gas. Further 
information on the performance of the tracking and dE/dx measurements can be found in 



reference |jTI|. These detectors are located inside a solenoid providing a magnetic field of 
0.435 T. Outside the solenoid are a time-of-flight scintillator array and a lead-glass electro- 
magnetic calorimeter with a presampler. Including the endcap electromagnetic calorimeters, 
the lead-glass blocks cover the range | cos6'| < 0.98. The next layer is the hadron calorimeter, 
consisting of the instrumented return yoke of the magnet. Finally, the detector is covered 
by several layers of muon chambers. In total, at least seven, and in most regions eight, 
absorption lengths of material lie between the interaction point and the muon detectors. 
This material is sufficient to absorb low-momentum muons produced at the vertex, but most 
muons with momenta above 2 GeV/c reach the muon detectors. The muon chambers are 
constructed as two different detector subsystems in the barrel and endcap; together, they 
cover 93 % of the full solid angle. 

This analysis uses events recorded at centre-of-mass energies within 3 GeV of the Z° peak 
during the 1992 — 1995 running period when the silicon microvertex detector was fully oper- 
ational. A total of 7.15 million multihadronic Monte Carlo events, and 4.88 and 1.93 million 
simulated bb and cc hadronic decays are also analysed. All samples were produced with 
the JETSET 7.4 Monte Carlo generator ||12|, with the fragmentation function of Peterson 
et al. |l^ for heavy quarks. The ACCMM model IQ tuned to the CLEO data [jl5| is used 



to describe the lepton momentum spectrum in b — > £ and b — i> c -^ £ decays, as described 
in [|l^. The Monte Carlo parameters were tuned to describe the OPAL data [jl^. All 



simulated events were passed through the full OPAL detector simulation package ||T8 



3 Event selection and analysis method 



Standard hadronic event selection [19| and detector performance requirements are applied to 
select a sample of 3.35 million events where the primary vertex can be reconstructed. Each 
event is divided into two hemispheres by the plane perpendicular to the thrust axis and 
containing the interaction point. The thrust value of the event is required to be greater than 
0.8 to suppress contributions from events containing more than two jets, for example from 
hard gluon radiation or events with gluon splitting into a heavy quark pair. The polar angle 
of the thrust axis, 6'th, must satisfy |cos6'th| < 0.75, to ensure that the event is contained 
within the central barrel region. A total of 2.15 million events satisfy these requirements. 

Lifetime tagging techniques are used to suppress the contributions from non-bb events. 
Hemispheres are tagged as containing b hadrons ("b-tagged") using a neural network al- 



gorithm IQ. A cut is apphed to the network output, selecting a sample of b- hemispheres 
with a purity of 92% and an efficiency of 30%. The b purity for the hemisphere b-tag is 
extracted directly from the data, as detailed in the next section. 

A search for lepton candidates is made in the hemisphere opposite a b-tagged hemisphere in 



^The coordinate system is defined such that the 2:-axis follows the electron beam direction and the x-axis 
points towards the centre of the LEP ring. The polar angle is defined relative to the +z-axis, and the 
azimuthal angle 4> is defined relative to the +a;-axis. 



events containing one or two such hemispheres. By using leptons found in the hemisphere 
opposite the b-tagged hemisphere rather than in the b-tagged hemisphere itself, a sample 
is obtained which does not bias the relative fraction of the different b hadron species. In 
addition, this method avoids introducing significant correlations between b flavour tagging 
and lepton selection. 

Jets are formed from charged tracks and electromagnetic energy clusters unassociated with 
tracks using a cone algorithm ||2I| , with a minimum energy of 5 GeV and a cone radius of 550 
mrad. Jet shape and momentum variables are used in two neural networks, NNb^ and NNbc^, 
trained respectively to distinguish direct decays, b — > ^, and cascade decays, b ^ c ^ £, from 
all other lepton sources, collectively termed as backgrounds. Separate neural networks are 
trained for electrons and muons. Details of the training of these neural networks are given 
in Section |4.1|. The distributions of the neural network outputs are compared for the data 



and the Monte Carlo to determine the fractions of events from b — ^ £ and b — *> c — ;> £ decays. 
Contributions from direct c —^ I decays are suppressed by the b-tagging requirement. 

To determine the fraction of leptons coming from b — >■ £ and b — > c — > £ decays, a binned log- 
likelihood fit is performed which uses the shapes of the distributions of the neural network 
output variables. The fit also yields {xe)- The number of b — > £ and b ^ c — > ^ decays is ob- 
tained by multiplying the fitted fractions by the number of selected leptons, corrected by the 
lepton detection efficiencies, derived from Monte Carlo. Dividing the numbers of b — *^ £ and 
b ^ c — > ^ candidates by the number of b-fiavoured hemispheres yields the branching ratios 
BR(b — ^ £X) and BR(b — *> c ^ DC) . The number of b-fiavoured hemispheres is extracted 
from the data using a double tagging technique. 



3.1 Purity of the sample of b-tagged hemispheres 

The purity Pb of the b-tagged sample of hemispheres is extracted from the data using 
a double-tagging technique, thereby minimising systematic uncertainties. The number of 
hemispheres N^ passing the b-fiavour tagging criteria described above is counted together 
with the number of events A^tt where both hemispheres are b-tagged. With the b-tagging 
efficiencies for the b, c and light flavours given by r/b, ?7c and r^uds, one can write 

Nt = 2NMu[RhVh + RcVc + (1 - i?b - Rc)Vnds)], (2) 

iVtt = iVMH [C^Rbvl + CcRcVc + Cuds ( 1 - Pb - Rc)vL] , (3) 

where A^mh represents the number of events that passed the multihadronic selection and 
flducial cuts. Pb and Pc are the fractions of multihadronic Z^ events decaying into bb 
and cc pairs, respectively. The hemisphere correlation coefficients Cq, where q represents 
any primary quark flavour, are given by the ratio Cq = Tj^q/r]'^, where rjq is the efficiency 
for tagging a hemisphere containing flavour q and rj^q is the efficiency for tagging both 
hemispheres. Equations ^ and |^ can be re-expressed in terms of the purity Pb instead of r/b 
using the definition Pb = 2A^MH-Rb^b/^t5 and then used to directly measure the b purity. 
Whilst the b purity can be determined from the direct solution of either Equation ^ or 
Equation |^ separately, the value of Pb is extracted by maximising the log-likelihood of both 
equations simultaneously to obtain the maximum statistical sensitivity. 



Input parameter 


value 


5Pb 


Re 
Rb 

Vc 

Vnds 
Cb 


0.177 ±0.008 

0.2169 ±0.0012 

0.0209 ± 0.0002 (MC stat.) ± 0.0014 (syst.) 

0.0034 ± 0.0000 (MC stat.) ± 0.0003 (syst.) 

1.0493 ± 0.0052 (MC stat.) ± 0.0052 (syst.) 


±0.19% 
±0.01% 
±0.34% 
±0.22% 
±0.01% 


total 




±0.45 % 



Table 1: Contributions to the systematic uncertainty 6Ph on the h purity Ph- For 
the input parameters taken from the Monte Carlo, the statistical and systematic 
uncertainties are shown separately. An error of 0.0000 indicates a value less than 
0.00005. 



A total of A^MH = 2 150 423 events passed the multihadronic event selection. These events 
contained At = 303 366 b-tagged hemispheres, of which Att = 45 351 have also the other 
hemisphere b-tagged. The world average values of i?b and Re are used as inputs. To 
extract Pb from Equations |^ and ^, the charm and light-flavour efficiencies 7]^ and r^u^s and 
the correlation for b events, Cb, are taken from Monte Carlo, while Cc and Cuds? which have 
negligible impact on the b purity determination, are taken to be unity. From the fit to the 
data, the purity of all b-tagged hemispheres is measured to be (91.901 ± 0.016) %, where the 
error is statistical. 

Extensive studies have been presented in a previous OPAL analysis on the systematic differ- 
ences between data and Monte Carlo for Cb, r/c and ?7uds [l0|. Similar studies were conducted 
for this analysis. This was performed separately for the data collected in 1992 (19% of the 
total data sample) where only r — information was available from the silicon microvertex 
detector, and for the larger part of the data, where z information was also available. The 
corresponding weighted averages for the combined sample are used as input parameters to 
the purity ffi. The systematic uncertainties on Pb are summarised in Table ^ The purity is 



Pb = (91.90 ± 0.02 (stat.) ± 0.45 (syst.)) % 



(4) 



for a b-tagging efficiency r^b of around 30 



3.2 Electron identification 



High-momentum electrons are searched for in the hemisphere opposite to a b-tagged hemi- 
sphere. Electrons are identified using an artificial neural network pOl. The six inputs 



used by this neural network are: the momentum and polar angle of the track; the ra- 
tio E/p of the electromagnetic energy and track momentum; the number of electromag- 
netic calorimeter blocks contributing to the energy measurement; the normalised ioniza- 
tion energy loss di?/da;|norm and its error. The normalised dE/dx value is defined as 
di?/da;|norm = {dE / dx — dE / dx\o) / ao) , where dii^/dx|o is the ionization energy loss expected 
for an electron and Uo the error. 



Candidate tracks must have a minimuin of 40 jet chamber hits usable for the determination 
of the energy loss, out of a possible 159 hits. Electrons are required to have momentum 
greater than 2 GeV/c. The momentum cut is applied to reduce the fraction of fake electrons 
and to restrict the analysis to the region where the input variables used for this neural 
network are properly modelled. Electrons from photon conversion and electrons from Dalitz 



decays are rejected as in [gO 



In the data sample, 301303 b-tagged hemispheres were selected when the sub-detectors re- 
quired for electron identification were fully operational and 29 516 electrons were found in 
the opposite hemispheres. Monte Carlo events are used to determine the electron identifica- 
tion efficiency after b-tagging for b — *> e and b — i> c — ;> e decays. These are measured to be 
eb^e = 0.5662 ± 0.0012 and eb^c^e = 0.3306 ± 0.0012, where the errors come from Monte 
Carlo statistics. The effect of the momentum cut at 2GeV/c is taken into account. 



3.3 Muon identification 

All hemispheres opposite a b-tagged hemisphere are searched for muon candidates. Muons 
are identified by associating central detector tracks with track segments in the muon detectors 



in two orthogonal planes [^. The muon candidates are also required to have momenta 
greater than 2GeV/c. In total, 44832 muons are found in the hemispheres opposite 302 577 
b-tagged hemispheres selected when the muon chambers were functional. The efficiencies for 
identifying muons in b ^ /i and b ^ c ^ yU decays with momentum greater than 2 GeV/c 
are eb-+^ = 0.6794±0.0011 and Cb^c^^ = 0.4277±0.0013, where the errors come from Monte 
Carlo statistics. 



4 Composition of the lepton sample 

The main contributions to the sample after selecting a lepton in a b-tagged event come from 
direct decays, b — > £, and cascade decays, h —>■ c ^ i. All remaining sources are collectively 
referred to as background. 

Fake muons form the largest source of background in the muon sample. These fake muons 
are mostly due to light mesons passing through the hadronic calorimeter without showering. 
Fake electrons are less common and consist mostly of mis-identified pions. These fake leptons 
tend to have less transverse momentum than leptons from either b— s>£orb— >c— >£ decays. 
True leptons in bb events not originating from the semileptonic decay of a b or c quark, for 
example electrons from photon conversions, are grouped together as non-prompt leptons. 
The decays h ^ c —>■ i, where the c comes from the virtual W boson decay, form an important 
background to b ^ c — > £ decays. The lepton coming from either of the two cascade decays 
tends to be produced with less transverse momentum with respect to the jet axis than a 
lepton coming from a direct h —>■ i decay. The selected lepton samples also contain small 
contributions from b — >■ r ^ £ and b — > J/ip — ^ i^i~ decays. A smaller contribution to the 
backgrounds comes from leptons from charm and light fiavour decays. Due to the high b 
purity of the selected data sample, this source is greatly suppressed. Finally, a very small 



Monte Carlo lepton candidates 


104 653 e 


173 735 /i 


b^£ 


54.7% 


43.0% 


b^c^£ 


27.1% 


23.0 % 


fake leptons 


3.3% 


18.5 % 


non-prompt leptons 


5.3% 


7.3% 


b^c^£ 


3.5% 


3.1% 


h^T^E 


2.5% 


1.8% 


b ^ J/^ -> i+t 


0.9% 


0.7% 


primary cc events 


2.4% 


2.2% 


primary uds events 


0.2% 


0.2% 


g^cc 


0.2% 


0.1% 


g^bb 


< 0.001 % 


< 0.001 % 



Table 2: Composition of the Monte Carlo sample of electron and muon candidates 
opposite a b-tagged hemisphere showing the contributions from b— >£, b— ^^c^^ 
and background. 



contribution comes from gluon splitting into cc and bb. The sources of background found in 
the Monte Carlo sample are shown in Table |^. 

Instead of attempting to reject these backgrounds, a fit for the fractions of b — *> £ and 
b ^ c — s> £ decays in the sample is performed using the two-dimensional distribution of the 
output variables of two neural networks: The first neural network, NNb^, is trained to sep- 
arate b — ^ ^ events from all other events while the second, NNbc£, is trained to separate 
b — i> c — *> £ events from all other events. Each of the neural networks is trained separately 
for electron and muon samples since the background is different in the two. 



4.1 Neural network training 

Jet and lepton kinematic variables in b-tagged events are used as input variables to train the 
neural networks NNb£ and NNbc^ to select direct and cascade leptons coming from b decays. 
By combining the information from several variables, more discrimination power is obtained 
than by using momentum information alone. About 90 000 muon and 70 000 electron candi- 
date tracks reconstructed in Monte Carlo events opposite a b-tagged hemisphere, were used 
to train NNb^ and NNb^- Eight kinematic variables were used, which are shown in Figure 



lepton momentum; 

lepton pt: the lepton transverse momentum calculated with respect to the nearest jet 
axis, excluding the lepton candidate itself; 

lepton jet energy: the energy of the jet containing the lepton; 



• sub-jet energy: the energy of the sub-jet (defined below) containing the lepton; 

• pt sum: the scalar sum of transverse momenta of all charged tracks in the lepton jet; 

• impact parameter: the impact parameter significance of the candidate lepton track 
with respect to the primary vertex; 

• lepton Qjeu the lepton charge multiplied by the jet charge (defined below) of the jet 
containing the lepton, including the lepton candidate track; 

• opposite Qjet'- the lepton charge multiplied by the jet charge of the most energetic jet 
in the hemisphere opposite the lepton. 

In h —>■ i decays, the lepton momentum spectrum refiects the hard fragmentation of the 
primary b hadron and is thus efficient at separating these leptons from other sources. The 
boost along the b jet direction results in a higher lepton momentum for b ^ £ than for other 
decays. Similarly, the decaying b hadron imparts more pt to the lepton in b — > ^ decays 
than in b — >• c — > £ decays. 

The total energy of the jet has sensitivity to leptons from direct decays since b jets are 
expected to have lower visible energy in semileptonic decays due to the emission of an 
energetic neutrino. 

The smaller mass of c hadrons relative to b hadrons forces the non-leptonic decay products 
from a charm decay to follow the lepton direction more closely. The neutrino in a charm 
decay also carries less energy on average than the neutrino in a primary b — *> £ decay. These 
differences mean that the energy deposited by neutral and charged particles in the vicinity 
of the lepton candidate, the lepton sub-jet energy, will be lower in b — *> £ decays than in 
b — i> c ^ £ decays. The lepton jet is therefore divided into two sub-jets |^, where the 



initial sub-jet seeds are the lepton track and the other tracks in the jet. Each track and 
unassociated electromagnetic cluster is then reassigned iteratively until each one is closer in 
angle to its assigned sub-jet axis than to the other. No track or cluster is added to the sub- 
jet containing the lepton beyond an invariant mass upper limit of 2.5GeV/c^. The "sub-jet 
energy" used for the neural network input refers to the sub-jet including the lepton. 

The scalar sum of pt for the jet characterises both the angular width and the multiplicity of 
the jet, both of which are known to differ significantly between b and charm and light-quark 

jets pi. 



The lepton impact parameter significance is the distance of closest approach of the track to 
the primary vertex divided by the uncertainty on this distance. Larger impact parameter 
significances are expected for leptons from decays such as b — ;> £ and h -^ c -^ i decays, than 
for tracks from the primary vertex such as fragmentation tracks. 

The final two variables consist of the reconstructed lepton charge multiplied by the jet 
charge for the jet associated with the lepton and for the most energetic jet in the b-tagged 
hemisphere. The jet charge is defined as 

^iei = ^^ o^g (5) 

2^P^ 



where Qi is the track charge, Pi is the track momentum and the sum runs over the charged 
tracks in the jet including the lepton candidate. Leptons from b — *> £ decays have the same 
charge as the weakly decaying b quark and thus the lepton Qjet variable shows a positive 
correlation between the lepton charge and associated jet charge. Leptons from b — »• c — ^ £ 
decays have opposite charge to the decaying b quark and hence show a negative correlation 
with the lepton jet charge. Leptons from b ^ c — > £ decays have a positive correlation with 
the lepton jet charge. In the absence of B" — B" mixing, the correlations between the lepton 
charge in one hemisphere and the jet charge in the opposite hemisphere, embodied in the 
opposite Qjet variable, are opposite to those of the jet associated with the lepton. 

The distributions of the eight input variables used by the neural networks are shown for 
electrons in Figure |^. The same good agreement between data and Monte Carlo simulation 
is found for muons for all input variables. Combining the information from these variables 
using neural networks allows not only to increase the separation power but also to include 
the correlations between the input variables. 

The distribution of the neural network outputs from NNb^ and NNbc^ for muons are shown 
in Figures |^ and ^. The results are shown for the three categories of Monte Carlo events: 
direct b decays b — >• £, cascade decays b ^ c ^ £, and the background. As an indicator of 
the goodness-of-fit, the x^/bin is 1.00 in Figure |^ and 0.88 in Figure ^ and is calculated 
using the statistical, systematic and modelling errors (see Section P). 

As can be seen from Figure ^, the neural network NNb^ gives good separation of b — > £ decays 
from all other types of decays found in the b-tagged sample; there is much less separation 
power between b ^ c ^ £ events and the background. The second neural network, NNbc^, 
gives further discrimination between b ^ c — > £ events and all other events as can be seen 
in Figure H. 



5 Calculation of the semileptonic branching ratios 

The semileptonic branching fraction is given by 

jVb Cb^^ jVb-tags ■ -T b 

where A^b-^^ is the number of hemispheres containing a semileptonic b decay and A^b the total 
number of true b hemispheres. The fraction of b ^ £ events determined by the fit, /(b — * £), 
times the number of lepton candidates, A^, yields the number of b — ^ £ decays in the selected 
lepton sample. The total number of b events that decayed semileptonically, N]^^i, is then 
obtained by correcting this number to account for the lepton detection efficiency, e^^e. The 
number of true b decays in the b-tagged sample, A'b, is obtained from the total number 
of b-tagged hemispheres, A'b-tags scaled by the sample purity, Pb, as extracted from the 
data. This calculation is performed separately for electrons and muons. By replacing all 
occurrences of b — > £ by b — > c — >■ £ in Equation ^, one obtains the corresponding equation 
for BR(b ^c^iX). 
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Figure 1: The eight input variables used by the neural networks NNb^ and NNbc^ 
to separate b — *> £ and b ^ c — *> £ decays. The contributions from the three classes 
of Monte Carlo events; h^i,h^c^£ and backgrounds, are superimposed and 
compared to the data after normalisation and rescaling using the fractions f{h—*£) 
and /(b ^ c — > 
only. 



derived in Section \5.2. The results are shown here for electrons 
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Figure 2: Distribution of the output variable for the neural network NNb^ trained 
specifically to distinguish b — *> £ events from all other types of events. The three 
categories b -^ £, b ^ c ^ £, and background are described in the text. The results 
are shown here for niuons only. 
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Figure 3: Distribution of the output variable for the neural network NNb^ trained 
specifically to distinguish b — *> c — *> £ events from all other types of events. The three 
categories b — * £, b ^ c ^ £, and background are described in the text. The results 
are shown here for niuons only. 
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5.1 Determining the lepton sample composition 

The fractions oi h —>■ i and h ^ c —>■ i decays are extracted from the data by performing a 
binned log-hkehhood fit to the data distributions of the neural network outputs using the 
shapes for each contribution obtained from the Monte Carlo. Two-dimensional distributions 
are formed using the neural network outputs NNb^ and NNbcf with 20 by 20 bins. There are 
five fit parameters: the Peterson b fragmentation model parameter e?,, /(b —>■ e), /(b ^ /i), 
/(b ^ c ^ e) and /(b ^ c ^ /i). The remaining fraction of selected candidates is assigned 
to the background sample: h —>■ c ^ i leptons, non-prompt leptons and fake leptons in bb 
events, and leptons found in charm and light flavour events. Their fraction is fixed to be 
1 — /(b ^ £) — /(b — i> c — i> £) in the fit for electrons and muons separately. The Peterson b 
fragmentation model parameter e^ is allowed to float in the fit, with a common value used for 
both electrons and muons, to greatly reduce the contributions to the systematic uncertainties 
from the modelling of this parameter. This is done by reweighting the Monte Carlo events 
such as to produce the desired fragmentation function. 



5.2 Results 

In total, 29 516 electrons and 44832 muons are selected in the data opposite to a b-tagged 
hemisphere. The results of the fit are /(b — > e) = 0.5726 ± 0.0042 and /(b ^ c ^ e) = 
0.2596 ±0.0055 for electrons and /(b ^ /i) = 0.4620 ±0.0034 and /(b ^ c ^ /x) = 0.2166 ± 
0.0051 for muons where the uncertainties are statistical. The Peterson fragmentation pa- 
rameter is determined to be eb = 0.00573 ± 0.00040(stat.). These results include some small 
corrections described in the next section. Using these corrected parameters, the fit can ade- 
quately reproduce the observed distributions of the neural network output variables in data 
for all selected leptons, as shown above in Figures ^ and ^. 

The correlation coefficient between /(b —>■ €) and /(b ^ c ^ £) is —0.35 for electrons and 
—0.26 for muons. The correlations of /(b — ^ e) and /(b -^ fi) with the common fragmenta- 
tion parameter are 0.32 and 0.39 respectively, whilst the fragmentation parameter produces 
a 0.13 correlation between /(b — > e) and /(b — > /i). The size of all other correlations is below 
0.06. 



From Equation y using the b purity Pb as determined in section |3.1| , the semileptonic branch- 
ing fractions are then determined to be 



BR(b- 


-^eX) = 


= (10.78 ±0.08)% 


BR(b- 


.fiX) - 


= (10.96 ±0.08)% 


BR(b -^ c - 


-^eX) = 


= (8.37 ±0.18)%, 


BR(b ^ c - 


./iX) = 


= (8.17 ±0.19)%, 



where the errors are statistical. 
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6 Systematic and modelling uncertainties 

In this section, possible systematic shifts and uncertainties on the five fit parameters are 
estimated by studying detector modelhng and experimental systematic effects on the results. 
The corresponding systematic uncertainties on BR(b — > iX), BR(b — > c — > ^X) and (xe) are 
compiled in Table |^. 



6.1 Modelling the semileptonic decay lepton spectrum 

Previous BR(b — *> iX) measurements @-0 depended heavily on the proper modelling of the 
semileptonic decay spectrum. The exact shape of the lepton momentum spectrum is not 
known and little theoretical progress has been made in recent years. The use of neural net- 
works to separate b — >■ £ and b — >■ c — > £ events from the background reduces the dependence 
of the branching ratio measurements on the shape of the lepton spectrum by making use of 
extra information such as jet shape and charge correlations. Nevertheless, the simulation of 
the b hadron decays and the prediction of the lepton momentum spectrum is still a large 
source of uncertainty. 

Different decay models are used to estimate the effects of the modelling on the fitted fractions. 
The events are reweighted to reproduce the lepton momentum spectrum in the rest frame 
of the b hadron as predicted by the different models. The ACCMM model [0] predictions 



are used for b — *> £ and b — i> c ^ £ decays to calculate the central values of BR(b — *> iX) 
and BR(b ^c^iX), using the prescription in |T§. The ISGW ^ and ISGW** modelsQ 



provide the iblcr deviations for the h —* i decays [|T^. Although these models were derived 
using only B" and B^ mesons, all b hadrons are reweighted. This has a very small effect on 
the central value but gives a more conservative estimate of the uncertainty on the modelling 
error than when only the B° and B^ decays are reweighted. The agreement between the 
data and these and other semileptonic decay models is further investigated in the Appendix. 

For the cascade decays, h -^ c ^ i, we use variations based on the CLEO measurements of 
the b — *> D spectrum combined with the ACCMM prediction for c — *> £ decays as described 
in [jl6[. Since the b-tagging requirement highly suppresses contributions from charm and 
light flavours, the error from the semileptonic decay modelling in charm events is negligible, 
and these are simply reweighted to the central ACCMM model as described in |16|]. Helicity 



effects in D* decays are not taken into account but are expected to be very small. 

Because of the minimum momentum cut of 2 GeV/c imposed on the selected leptons, the 
measured lepton identiflcation efficiencies correspond to a restricted momentum range. The 
effect of the extrapolation below the minimum momentum cut-off is taken into account when 
evaluating the lepton selection efficiency corresponding to the different models. 



^This corresponds to a modification of the ISGW model introduced by the CLEO collaboration whereby 
all P-wave contributions (the D** contributions) are increased from the nominal 11 % derived in the original 
model to 32 % to better describe their data. 
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6.2 Pragmentation models 

Several models have been proposed to describe the heavy quark fragmentation process. The 
function of Peterson et al. |13| was used to simulate fragmentation in bb and cc events 
in the Monte Carlo. For b hadrons, the Peterson parameter is determined from the fit by 
reweighting the Monte Carlo events. For c hadrons, the parameter is varied between 0.025 to 
0.031 to obtain a {xe) for charm hadrons of 0.484±0.008 ||T6[. As suggested in |T^, the models 



of Collins and Spiller ||2^, and Kartvelishvili et al. [^ are used to estimate the systematic 
uncertainties coming from the shape of the b quark fragmentation function, quoted as the 
+ lcr and —la errors respectively. These models also have one free parameter. The Monte 
Carlo is reweighted to simulate each function and the free parameter as determined from 
the fit. The effects of the differing fragmentation functions and fitted parameters on the 
lepton efficiencies are also included in the fragmentation modelling error. The systematic 
uncertainties associated with the b fragmentation models are determined from the observed 
variations in the derived values of branching fractions obtained with the various functions. 
The agreement between the data and these functions is further discussed in the Appendix. 



6.3 Lepton detection efficiencies 

The electron detection efficiencies in data and Monte Carlo are compared for two pure 
electron samples: a low-momentum electron sample coming from photon conversions, and 
a high-momentum electron sample from Bhabha events. From these comparisons and from 
studying the neural network input variable distributions, a systematic uncertainty on the 
electron detection efficiency of 4%, and an uncertainty on photon conversion rejection of 
0.8 % is obtained 



The muon detection efficiency in data and Monte Carlo are compared for two samples: muon 
pairs from two-photon interactions and muons from Z° -^ f^^fJ'^ decays. The first sample 
yields muons in the momentum range of 2 to 6 GeV/c, while the second sample gives muons 
with momentum greater than 30 GeV/c. In the data, 57 % of the selected muon sample is 
found in the lower range while only 0.2 % has momentum above 30 GeV/c. By comparing the 
selected muon samples, the Monte Carlo was found to slightly underestimate the efficiency 
in the data and so a multiplicative correction factor of 1.013 is applied to the measured 
efficiencies. A systematic uncertainty of 1.9 % is assigned to the muon detection efficiencies. 

The systematic uncertainty associated with the momentum extrapolation below 2 GeV/c is 



discussed in Section 3.1. 



6.4 Detector resolution effects 

The tracking resolution and reconstruction efficiency could be slightly different for data and 
Monte Carlo. The reconstructed track parameters are smeared by ±10 % in the Monte Carlo 
and the lepton detection efficiencies and the fit fractions recalculated. The difference from 
the final values for BR(b —>■ iX) , BR(b ^ c ^ iX) and the b fragmentation parameter are 
used as estimate of this source of systematic uncertainty. 
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6.5 Fake lepton rates 

To study the fake rate in the muon sample, three different samples are used: identified pions 
in Kg — * TT^TT" decays, three prong r decays, and a kaon enriched sample based on dE/dx 
requirements. From a comparison of the fake muon rates in data and Monte Carlo for these 
samples, it was determined that a correction factor of 1.11 it 0.12 must be applied to the 
Monte Carlo events. Accordingly, the relative fraction of fake muon events is changed by 
+11 % using a reweighting technique. The weights are varied by ±12 % and the data refitted. 
For electrons, studies on the fake rates were conducted using samples of photon conversion 



electrons and Bhabha events pO[. The fake rates measured in the data and Monte Carlo were 
consistent, such that no correction was required, but an uncertainty of ±21 % is assigned for 
the fake electron rate. 



6.6 b tagging purity 

The systematic uncertainty on the value of the b purity obtained from the data is discussed in 
Section p.l| . This constitutes a 0.49 % relative error on the final values for BR(b — ^ £X) and 
BR(b —>■ c ^ iX). The errors on the b purity given in Table |I]have been further subdivided 
in Table |^ to show separately the contributions from uncertainties in Re, Rh, {xe)c (the 
mean fraction of the beam energy carried by the weakly decaying charmed hadrons), gluon 
splitting to bb and cc pairs, the branching fraction of charmed mesons into K^, charmed 
mesons lifetimes, decay multiplicities of charmed mesons and charm production fractions. 
The errors resulting from the uncertainty in the b purity due to detector resolution and finite 
Monte Carlo statistics are combined with the other contributions from these sources of error. 



6.7 Other sources of systematic uncertainties 

Several other sources of systematic uncertainties have been investigated. The Monte Carlo 
is reweighted to simulate the desired parameters and the fit is repeated to assess the contri- 
butions to the systematic uncertainty. These sources are summarised in Table |^. Since their 
effects on the branching ratios and (xe) are small, they are only described briefly here. 

Finite Monte Carlo sample size: This includes contributions from the evaluation of rjc, 
r^uds and Cb from Table |l]. 

b hadron species: The production fraction for Ab baryonsQ is set to (10.1 1ll'i) % [^ and 
the semileptonic branching fraction for inclusive Ab is taken to be (7.4 ± 1.1) % P,|28[|. 

b — *■ Xuii^ transitions: The b -^ Xuii^ branching fraction is set to (1.84 ± 0.79) x 10~^, 



the combined value of ||29| 



B° — B° mixing: The B^ — B[J mixing parameter is set to Xd = 0.172 ± 0.010 [^ whilst 
50% Bg — B° mixing is assumed [Q. 

■^Ab denotes all weakly-decaying b baryons produced in Z" decays 
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Fake lepton spectrum: Small variations in the distributions of fake leptons are allowed 
by shifting the momentum spectrum of fake leptons and non-prompt lepton in bb events by 
25MeV/c, roughly ±0.5% of the mean momentum for fake and non-prompt leptons in bb 
events. 

Contributions from b — > t — > £: The h ^ t ^ i branching fractions are set to (0.463 ± 
0.071) % and (0.452 ± 0.069) % 0] for electrons and muons, respectively. 

Contributions from b — > c — > i: This branching fraction is set to (1.62 ILq'^q) % as 
derived in [ pO| . 

Contributions from b — > J/V' — ^ £'^£~: The experimental value given in for BR(b -^ 

J/ip) is (1.16 ± 0.10) %. Combined with a recent BES measurement for BR.{J/iIj — > i'^i") = 
(5.87 ± 0.10) % 0, this gives BR(b -^ J/^ -^ £+£") = (0.0681 ± 0.0060) %. 

Effect of Ab polarisation: Leptons coming from semileptonic Ab decays are reweighted to 
simulate a spectrum corresponding to —56% polarisation according to ||32|. The systematic 
errors are calculated using the polarisation range —13% to —87%, the 95% confidence level 
limits from |3^ . 



6.8 Consistency checks 

To test the fitting procedure, the Monte Carlo sample is divided into two equal sub-samples. 
The first sample is used as a substitute for the real data in the fit while the second is 
retained as the Monte Carlo sample. The fitted fractions for the first sample can then 
be compared to the true information from the Monte Carlo. The fitted fractions and the 
Peterson fragmentation parameter all agree with the true values within statistics. 

Various tests are performed on the data to check the stability of the results by varying the 
selection criteria. Firstly, the effect of changing the minimum lepton momentum requirement 
on the measured value for BR(b — *> iX.) is investigated. The minimum lepton momentum cut 
is increased from the nominal 2.0 GeV/c to 5.0 GeV/c in steps of 0.5 GeV/c and BR(b — ^ £X) 
is recalculated each time. This test was found to yield good agreement when performed 
using the Monte Carlo test samples discussed above. Similarly, a cut is applied on the 
neural network NNb£ output variable shown in Figure |[ The cut is increased by steps of 0.1 
from 0.0 to 0.6. A similar test is performed to check the stability of the BR(b — > c — >^ iX.) 
results by imposing a cut on the neural network NNbc£ output variable shown in Figure ^ 
The cut is increased up to a value of 0.4. For all these tests, the variations observed are 
found to be statistically consistent with the central values calculated for BR(b — > £X) and 
BR(b^c^£X). 

Varying the binning used to perform the fit has no significant influence on the central values 
for BR(b — > £X) and BR(b — > c -^ £X). The central results are derived using 20 by 20 bins 
for the 2-dimensional fit to NNb£ and NNbc^. This range is varied from 5 by 5 up to 40 by 
40 bins, yielding consistent values for the branching fractions. 

Lastly, the data are divided into four sub-samples corresponding to the different years in 
which the data were taken. The b-tagging purity is recalculated for each data sub-set 
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Parameter 


BR(b -^ e) 


BR(b ^ c ^ e) 


BR(b -^ fi) 


BR(b ^ c ^ ^) 


(xe) 


model-dependent sources 


b->£ 
fragmentation 


-0.078 
+0.207 
-0.072 
+0.057 
+0.047 
-0.028 


+0.126 
-0.211 
+0.149 
-0.059 
+0.225 
-0.144 


-0.101 
+0.221 
-0.064 
+0.058 
+0.096 
-0.070 


+0.206 
-0.320 
+0.168 
-0.048 
+0.236 
-0.180 


-0.0051 
+0.0081 
+0.0009 
-0.0008 
-0.0118 
+0.0102 


total models 


+0.220 
-0.110 


+0.298 
-0.262 


+0.248 
-0.139 


+0.355 
-0.370 


+0.0131 
-0.0129 


systematic sources 


lepton efficiencies 


TO.440 


^0.341 


TO.208 


TO.155 




detector 


±0.074 


±0.113 


±0.055 


±0.086 


±0.0004 


lepton fake rates 


±0.006 


TO.048 


±0.037 


TO.106 


T0.0003 


Ph ■ Re 


±0.022 


±0.017 


±0.022 


±0.017 




Rh 


to.ooi 


to.ooi 


TO.OOI 


TO.OOI 




{xe)c 


±0.004 


±0.003 


±0.004 


±0.003 




g^bb 


±0.016 


±0.013 


±0.016 


±0.013 




g^CC 


±0.010 


±0.008 


±0.010 


±0.008 




BR(D^ K^) 


±0.011 


±0.008 


±0.011 


±0.008 




D° lifetime 


±0.002 


±0.002 


±0.002 


±0.002 




D^*" lifetime 


±0.003 


±0.002 


±0.003 


±0.002 




Dg lifetime 


±0.001 


±0.001 


±0.001 


±0.001 




D charged mult. 


±0.011 


±0.008 


±0.011 


±0.008 




D neutral mult. 


TO.024 


TO.018 


TO.024 


TO.018 




/(c^D±, DO) 


±0.017 


±0.014 


±0.017 


±0.014 




/(c ^ Ds) 


±0.002 


±0.001 


±0.002 


±0.001 




/(c ^ Ac) 


TO.007 


TO.005 


TO.007 


TO.005 




MC statistics 


±0.019 


±0.042 


±0.022 


±0.049 


±0.0010 


b hadron species 


TO.013 


±0.022 


TO.012 


±0.030 


T0.0006 


b -^ Xuiv 


±0.004 




±0.009 


±0.022 


TO.0020 


B mixing 


±0.002 


±0.016 


TO.002 


±0.007 


±0.0002 


fake lepton spectrum 




^0.003 


TO.002 


TO.042 




h^T^e 


TO.026 


TO.013 


TO.021 


TO.019 


±0.0003 


h^c^i 


TO.004 


^0.081 


TO.023 


TO.064 


±0.0003 


j/V' ^ i^r 


TO.004 




TO.002 




±0.0001 


Ab polarisation 


±0.004 


±0.006 


±0.005 


±0.026 


-0.0013 
+0.0020 


experimental systematic 


±0.450 


±0.377 


±0.227 


±0.234 


+0.0031 
-0.0027 



Table 3: Summary of all model-dependent and experimental systematic uncertainties 
on BR(b -^ iX) and BR(b ^ c ^ iX) (shown separately for electrons and muons), 
and {xe)- All errors are absolute errors given in percent (except for {xe))- The sign 
on each contribution indicates the correlation between this systematic uncertainty 
and the final results. 
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separately. Taking the uncorrelated systematic errors into account, all sub-samples are 
found to be statistically consistent with the full data sample. 



7 Results and conclusions 



All the relevant quantities needed to calculate the branching ratios BR(b -^ iX) and 
BR(b ^ c — i> ^X) are given in Table |^. The values 

0.22 



BR(b- 


-^eX) 


BR(b- 


./iX) 


BR(b ^ c - 


-^eX) 


BR(b ^ c - 


^liX) 



(10.78 ± 0.08 (stat.) ± 0.45 (syst.) Ig;^^ (model)) % 



0.11 
0.25 







(10.96 ± 0.08 (stat.) ± 0.23 (syst.) Z^Z (model)) % 



0.14 
0.30 







(8.37 ±0.18 (stat.) ±0.38 (syst.) ™ (model)) % 







-- (8.17 ± 0.19 (stat.) ± 0.23 (syst.) ^Jg? (model)) % 

are obtained for the semileptonic branching ratios for electrons and muons, consistent with 
lepton universality. These four branching ratios are combined together to obtain 

BR(b ^ iX) = (10.83 ± 0.10 (stat.) ± 0.20 (syst.) toll (model)) % 
BR(b ^c^iX) = ( 8.40 ± 0.16 (stat.) ± 0.21 (syst.) tl'S) (model)) % 



0, 



taking into account the full covariance matrix as in |T^. The BR(b — >■ iX) measurement is 
the most precise to date at the Z° resonance whereas BR(b —>■ c ^ iX) is more precise than 
the current world average value of (7.8±0.6)% [0]. The value derived for BR(b ^ c — > iX) is 
outside the range given by BR(b -^ c — > eX) and BR(b — > c ^ /iX) due to large off-diagonal 
terms in the covariance matrix and strong correlations with the BR(b -^ iX) measurement. 
The statistical error on BR(b -^ iX) is larger than the individual errors on BR(b -^ eX) and 
BR(b — >■ /iX) since the statistical errors from BR(b — *> c — >■ eX) and BR(b — *> c ^ /iX) also 
contribute. The full systematic correlation matrix is given in Table ^ 

The BR(b — >■ eX) and BR(b -^ fiX) measurements presented here are consistent with the 
current average of all Z° measurements, BR^^ = (10.99 ± 0.23)% 0, based on a global 
fit to several electroweak parameters and including specific measurements of BR(b —>■ iX) 
0-^]. On the other hand, the measurement of BR(b — > iX) is still larger than the T(4S) 
measurement of BRf^;^ = (10.45 ± 0.21) % 0, when it is expected to be lower, as explained 
in Section |I]. If the lifetime ratio correction is applied, the difference between this result and 
the T(4S) measurement is about 1.8 standard deviations. 

This measurement is also consistent with theoretical calculations as can be seen in Figure ^ 
where a correction factor of 0.974 corresponding to the lifetime ratio Tb/re has been 
applied to both the theoretical calculations and the T(4S) value of BR(B — >■ Xi) to allow 
direct comparison with the Z^ results. No correction is applied to the values of (ric), the 
average number of charm hadrons produced per b decay. 

From the fitted fragmentation model parameters, the average value of the fraction of the 
beam energy carried by the weakly decaying b hadron is obtained, giving 

(xe) = 0.709 ± 0.003 (stat.) ± 0.003 (syst.) ± 0.013 (model) 
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electrons 


muons 


n 


0.9190 ± 0.0002 (stat.) ± 0.0045 (syst.) 


-^b-tags 


301303 
29516 


302577 
44832 


/(b ^ i) 


0.5662 ±0.0231 (syst.) 
0.5726 ± 0.0042 ± 0.0041 (syst.) 


0.6794 ±0.0129 (syst.) 
0.4620 ±0.0034 ±0.0031 (syst.) 


BR(b ^ £X) 


(10.780 ± 0.079 ± 0.450 t^J^) % 


(10.964 ± 0.081 ± 0.227 1[J:?^^) % 


/(b ^ c ^ £) 


0.3306 ±0.0135 (syst.) 
0.2596 ± 0.0055 ± 0.0047 (syst.) 


0.4277 ±0.0081 (syst.) 
0.2166 ± 0.0051 ± 0.0045 (syst.) 


BR(b ^ c ^ £X) 


(8.370 ± 0.177 ± 0.377 tl]:l^)% 


(8.167 ± 0.192 ± 0.234 ^^^f^^ % 



Table 4: Results for the data sample including all systematic uncertainties for elec- 
trons and muons. The uncertainties from detector resolution have been added to 
the errors on the fitted fractions /(b — > £) and f{h ^ c ^ i). The uncertainties due 
to semileptonic decay and fragmentation modelling are shown in the last error on 
the branching fractions. 





BR(b^eX) BR(b^c^eX) BR(b ^ ^X) BR(b ^ c ^ /iX) 


BR(b- 


-^eX) 


1.00 


BR(b- 


-^ c ^ eX) 


0.40 1.00 


BR(b- 


-.^X) 


0.34 -0.22 1.00 


BR(b- 


^ c ^ /iX) 


-0.26 0.53 -0.22 1.00 



Table 5: The full systematic correlation matrix from the averaging procedure of 
BR(b -^ eX), BR(b -^ fiX), BR(b -^ c ^ eXj and BR(b ^ c ^ fiX). 
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Figure 4: The number of charm hadrons per b hadron decay, {uc) , as a function of 
BR(b -^ iX) for the Z^ /|/ and BR(B -^ XE) for the T(4S) result /|/. The result 
derived in this analysis is superimposed as a vertical line with error bars. The shaded 
box represents the theoretical calculation in the framework of Heavy Quark Effective 
Theory including higher-order corrections as given in ^ for different assumptions 
for the ratio of the renormalisation scale fi and the b quark mass mi,, and the ratio 
of the c and b quark masses, m,c/m,f,. A correction ofry^/r^ has been applied to the 
theoretical prediction and to the T(4S) result, as described in the text. 
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where the modehing error is dominated by the choice of b fragmentation model. 

All the measurements presented here are statistically independent of, and consistent with 
similar results derived in a previous OPAL analysis [^ where the quantities BR(b — *> iX) = 
(10.5 ± 0.6 (stat.) ± 0.5 (syst.)) %, BR(b ^ c ^ £X) = (7.7 ± 0.4 ± 0.7) % and (xe) = 
0.697 ± 0.006 (stat.) ± 0.011 (syst.) were extracted from a fit for these and several other 
parameters (including i?b)- However, the uncertainties related to assessing the systematic 
correlations between these old results and those presented in this paper means that no 
overall gain in precision is obtained by combining them. Therefore the results presented 
here supersede the results previously published in 0. 
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Appendix 



Semileptonic b decay and fragmentation models 

For the first time at LEP, an attempt is made to probe the level of agreement between the 
data and various semileptonic decay models. This test is also sensitive enough to allow a 
closer examination of different fragmentation functions. In this section, several theoretical 
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h ^ i decay models are investigated. These models affect both the lepton total and trans- 
verse momentum spectra. For each h ^ i decay model, we use three different fragmentation 
functions, those of Peterson et al. [|13|, Collins and Spiller [^ and Kartvelishvili et al. [ p7 



These functions primarily affect the lepton total momentum spectrum, leaving the trans- 
verse momentum unchanged. The same models are used to simulate the cascade decays 
b ^ c ^ £, and to assess the modelling uncertainties, as described in Section [O] . 

The six b — > £ decay models investigated are: 



1. ACCMM model [|T4|: all parameters were tuned to the CLEO data |T5|. Their val- 



ues are fixed as given in [^: the Fermi momentum of the spectator quark, pp = 
298MeV/c, the mass of the charm quark, rric = 1673MeV/c^, and the mass of the 
spectator quark, mgp = 150 MeV/c 



2 



2. ISGW model [25|: this model has no free parameters and the D** contributions are 



predicted to account for 1 1 % of all b decays. 
3. ISGW** model: the ISGW model modified such as to allow the total contributions 



from D** to account for 32% of all b decays, as described in [|T6|. 



4. ISGW2 1^: a revised version of the ISGW model incorporating constraints from 
heavy quark symmetry, hyperfine distortions of wave functions, and form factors with 
more realistic high-recoil behaviour. This model predicts that the sum of all D** 
contributions accounts for 9.3 % of the total b decay width. The b -^ X„£z/ branching 
ratio was fixed to 2.2 % of all b decays and mixing was assumed between the t] and t] 
final states. 

5. ISGW2**: the ISGW2 model modified to allow the sum of all D** contributions to be 
a free parameter of the fit. Best agreement with OPAL data is found when the D** 
contribution amounts to (45 ± 3 (stat.) ± 3 (syst.) ± 4 (model)) % of the total width, 
instead of the 9.3 % set in the original model. The model error contains uncertainties 
from both the b fragmentation (following the same procedure as described in Section 
|6.2| ) and b ^ c ^ £ decay models. The Peterson fragmentation model is used to derive 
the central value. 

6. ACCMM*: the ACCMM model with free parameters. The Fermi momentum of the 
spectator quark, pp, and the mass of the charm quark, rric, are treated as free param- 
eters in the fit, giving pp = (837 ± 143 (stat.) ± 132 (syst.) tUt (model)) MeV/c and 
rric = (1287±100 (stat.) ±87 (syst.) till (model)) Me V/c^. The mass of the spectator 
quark is kept fixed at 150 MeV/c^ as in |T6[ . 

The semileptonic decay model parameters pp and rric strongly depend on the choice 
of the fragmentation model used. The Peterson model is used to derive the central 
values. The modelling errors given here again correspond to the b fragmentation and 
b ^ c — > £ decay model errors added in quadrature. The correlation coefficient between 
these model parameters is —0.970. The decay model parameters are consistent with 
the calculated pp value of about 550 MeV/c derived in [Q using the relativistic quark 
model, and the world average charm mass of 1100 to 1400MeV/c^ taken from Q. The 
b -^ Xudv branching ratio was fixed to 2.7% of all b decays. 
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Figures |^ and | show the fitted distributions for each of these models compared to the data 
in the b — ^ £ peak region of the NNb^ neural network output distributions for electrons and 
muons, respectively. The fit is performed over the full range of the neural network output 
(from zero to one). The region shown corresponds to NNb^ > 0.8, which, from the fit results 
to the data, is approximately 93 % pure in b — > £ decays. 

The results for BR(b —>■ iX) and BR(b — *■ c — *■ £X) obtained with these different models are 
summarised in Table ^ together with the statistical, systematic and modelling uncertainties. 
All errors are calculated according to the procedures outlined in the preceding sections, 
apart from the modelling error which accounts for b ^ c — > £ decay modelling only. The 
values obtained for the decay model parameters as well as for the free parameter in the 
fragmentation functions are also given in Table |^. The results for (xe) corresponding to the 
various fitted fragmentation functions are also listed. These can be compared to the value 
of (xe) = 0.702 ±0.008 from [Q obtained from a multi-parameter fit to several electroweak 
parameters. The x^/bin is calculated using the portion of the NNb^ output shown in Figures 
^ and P using statistical, systematic and modelling uncertainties from both the electron and 
muon samples. These are given as an indicator of the agreement between these models and 
the data. 

The accuracy of the test does not allow ruling out specific b — > £ models, although some 
trends are clear: 



The fragmentation functions of Peterson et al. and Kartvelishvili et al. provide equally 
good fits to the data. The fragmentation function of Collins and Spiller is generally 
disfavoured by our data. 

The semileptonic decay models ISGW** and ISGW2** best agree with our data when 
used with the b fragmentation models of Peterson et al. or Kartvelishvili et al. However, 
these models are less theoretically sound since modifications to the original models were 
needed to allow the overall fraction of D** contributions to be a free fit parameter when 
this fraction is in fact one of the predictions of the models. 

For the ACCMM* semileptonic decay model, the best fit to the data is obtained with 

Pf = (837 ± 143 (stat.) ± 132 (syst.) tfst (model)) MeV/c, 
me = (1287 ± 100 (stat.) ± 87 (syst.) tlH (model)) MeV/c^ 

when the mass of the spectator quark is kept fixed at 150MeV/c^. These results are 
derived using the b fragmentation model of Peterson et al. 

The ISGW2 model gives a worse agreement with our data than the ISGW model, with 
all fragmentation models. 
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Figure 5: The fitted distributions for the NNb^ neural network output for electrons 
with (a) the ACCMM, ISGW and ISGW** models; (b) the ISGW2, ISGW2** and 
ACCMM* models; (c) the ACCMM model with the Peterson, Collins and Spiller and 
Kartvelishvili fragmentation functions. The Peterson function is used to describe 
the fragmentation in (a) and (b). The shaded area shows contributions from sources 
other than h —* e in the data, as extracted from the fit. 
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Figure 6: The fitted distributions for the NNb^ neural network output for muons 
with (a) the ACCMM, ISGW and ISGW* models; (b) the ISGW2, ISGW2** and 
ACCMM* models; (c) the ACCMM model with the Peterson, Collins and Spiller and 
Kartvelishvili fragmentation functions. The Peterson function is used to describe 
the fragmentation in (a) and (b). The shaded area shows contributions from sources 
other than h —* fi in the data, as extracted from the fit. 
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model 


h ^^ £ model 
parameters 


Fragmentation 
parameter 


{^e) 


BR(b -^ eX) 


BR(b -^ fiX) 


BR(b ^ e ^ eX) 


BR(b ^ e ^ ^X) 


xVbin 


Peterson et al. 


ACCMM 


fixed 


0.00573 ± 0.00062 


0.709 ± 0.004 


10.78 ±0.08 ±0.45 ^[J;^^ 


10.96 ±0.08 ±0.23 


-0.06 
+0.06 


8.37 ±0.18 ±0.38 to,ll 


8.17 ±0.19 ±0.23 


+0.17 
-0.05 


64/48 


ISGW 


fixed 


0.00655 ± 0.00070 


0.705 ± 0.004 


10.70 ±0.08 ±0.45 ^°-°l 


10.86 ±0.08 ±0.23 


-0.07 
+0.06 


8.50 ±0.18 ±0.38 toll 


8.37 ±0.19 ±0.24 


+0.16 
-0.04 


98/48 


ISGW** 


fixed 


0.00456 ± 0.00051 


0.718 ± 0.004 


10.99 ± 0.08 ± 0.46 ^°-°l 


11.19 ±0.08 ±0.23 


-0.06 
+0.06 


8.16 ±0.18 ±0.37 talt 


7.85 ±0.20 ±0.23 


+0.19 
-0.07 


37/48 


ISGW2 


fixed 


0.00787 ± 0.00083 


0.698 ± 0.004 


10.69 ± 0.08 ± 0.45 ^f,':!]^ 


10.86 ±0.08 ±0.23 


-0.07 
+0.06 


8.62 ±0.18 ±0.38 torn 


8.55 ±0.19 ±0.24 


+0.15 
-0.03 


131/48 


ISGW2" 


/d,*=45±5 

PF= 837 tTr 
m, = 1287 t\tl 


0.00446 ± 0.00055 


0.719 ± 0.004 


10.95 ± 0.08 ± 0.45 ^"[^"^ 


11.15 ±0.08 ±0.23 


-0.08 
+0.07 


8.17 ±0.18 ±0.37 to.m 


7.85 ±0.20 ±0.23 


+0.40 
-0.08 


35/48 


ACCMM* 


0.00465 ± 0.00054 


0.717 ±0.004 


10.95 ± 0.09 ± 0.46 ^JJ;J^ 


11.15 ±0.09 ±0.23 


-0.10 
+0.08 


8.16 ±0.18 ±0.37 toM 


7.84 ±0.20 ±0.23 


+0.51 
-0.09 


38/48 


Collins and Spiller 


ACCMM 


fixed 


0.00342 ± 0.00062 


0.698 ± 0.004 


10.83 ± 0.08 ± 0.45 ^EJ;^^ 


11.06 ±0.08 ±0.23 


-0.07 
+0.06 


8.60 ±0.18 ±0.38 toll 


8.40 ±0.19 ±0.24 


+0.18 
-0.06 


148/48 


ISGW 


fixed 


0.00421 ± 0.00074 


0.693 ± 0.004 


10.74 ±0.08 ±0.45 '^[^ 


10.95 ±0.08 ±0.23 


-0.07 
+0.06 


8.72 ±0.18 ±0.38 toll 


8.61 ±0.19 ±0.24 


+0.17 
-0.05 


202/48 


ISGW** 


fixed 


0.00241 ± 0.00044 


0.705 ± 0.004 


11.05 ±0.08 ±0.46 -°-°l 


11.30 ±0.08 ±0.23 


-0.07 
+0.06 


8.39 ±0.18 ±0.38 tolt 


8.09 ±0.20 ±0.23 


+0.20 
-0.08 


84/48 


ISGW2 


fixed 


0.00556 ± 0.00096 


0.687 ±0.004 


10.72 ± 0.08 ± 0.45 ^^'^^f/g 


10.94 ± 0.08 ± 0.23 


-0.07 
+0.06 


8.83 ±0.18 ±0.39 toll 


8.79 ±0.19 ±0.24 


+0.16 
-0.04 


253/48 


ISGW2** 


/j3„ =43 ±5 

CTn +180 

PP = 679 Ii92 

™ _ 1987 +146 

nic - 1287 _j3g 


0.00251 ± 0.00044 


0.704 ± 0.004 


11.00 ±0.08 ±0.45 ^!J;[j? 


11.24 ±0.08 ±0.23 


-0.08 
+0.07 


8.43 ±0.18 ±0.38 torn 


8.14 ±0.20 ±0.23 


+0.41 
-0.06 


80/48 


ACCMM* 


0.00252 ± 0.00043 


0.704 ± 0.004 


10.94 ± 0.09 ± 0.45 ^"'H 


11.19 ±0.09 ±0.23 


-0.10 
+0.08 


8.43 ±0.18 ±0.38 toll 


8.14 ±0.20 ±0.23 


+0.49 
-0.09 


88/48 


Kartvelishvili et al. 


ACCMM 


fixed 


10.04 ±0.57 


0.720 ± 0.005 


10.75 ±0.08 ±0.45 ^[J'JJ^ 


10.89 ±0.08 ±0.23 


-0.07 
+0.06 


8.23 ±0.18 ±0.37 toll 


7.99 ±0.19 ±0.23 


+0.16 
-0.04 


41/48 


ISGW 


fixed 


9.40 ±0.54 


0.714 ±0.005 


10.69 ± 0.08 ± 0.45 ^°-°l 


10.80 ±0.08 ±0.23 


-0.07 
+0.06 


8.36 ±0.18 ±0.37 toll 


8.20 ±0.19 ±0.24 


+0.15 
-0.03 


56/48 


ISGW** 


fixed 


11.23 ±0.63 


0.729 ± 0.005 


10.94 ± 0.08 ± 0.45 ^EJ;^^ 


11.10 ±0.08 ±0.23 


-0.06 
+0.06 


8.01 ±0.18 ±0.37 tolt 


7.66 ±0.20 ±0.22 


+0.18 
-0.06 


48/48 


ISGW2 


fixed 


8.58 ±0.49 


0.706 ± 0.005 


10.69 ± 0.08 ± 0.45 ^°-°l 


10.81 ±0.08 ±0.23 


-0.07 
+0.06 


8.48 ±0.18 ±0.38 toll 


8.39 ±0.19 ±0.24 


+0.14 
-0.02 


73/48 


ISGW2** 


/d**=46±5 
PF = 1063 1^8 
nic = 1153 tilt 


11.44 ±0.67 


0.731 ±0.005 


10.91 ± 0.08 ± 0.45 ^J5;[j? 


11.07 ±0.08 ±0.23 


-0.07 
+0.06 


8.01 ±0.18 ±0.37 tolr 


7.65 ±0.20 ±0.22 


+0.39 
-0.09 


53/48 


ACCMM* 


10.96 ±0.64 


0.727 ±0.005 


10.95 ±0.09 ±0.46 ^EJ;J^ 


11.11 ±0.10 ±0.23 


-0.10 
+0.11 


8.00 ±0.18 ±0.37 toAO 


7.64 ±0.20 ±0.22 


+0.52 
-0.20 


39/48 



Table 6: Branching fractions BR(b — > ^X) and BR(b ^ c ^ iX.) (given in %) derived by comparing the data to three fragmentation 
functions and various semileptonic decay models for h ^ i decays. The uppermost hne corresponds to the central results, as 
discussed in section 7. The quoted errors on the branching ratios correspond to the statistical, systematic and b ^ c ^ .^ modelling 
errors, respectively. The fitted h ^> i decay model parameters are also given when appropriate. The fitted b fragmentation function 
parameters, and the corresponding values for (xe) are presented. The combined statistical, systematic and b — > c — > ^ modelling 
errors are given for all fitted model parameters. The x^/bin is calculated using the portion of the NNb^ output shown in Figures 5 
and 6, using all uncertainties from both the electron and muon samples; these are given as an indicator of the goodness-of-fit. All 
models and their parameters are described in the text 
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